Resolving Wave Propagation in Anisotropic Poroelastic Media Using Graphical Processing Units (GPUs)

Abstract Biot's equations describe the physics of hydromechanically coupled systems establishing widely recognized theory poroelasticity. This has a broad range applications in Earth and biological sciences as well engineering. The numerical solution is challenging because wave propagation fluid pressure diffusion processes occur simultaneously but feature very different characteristic time scales. Analogous to geophysical data acquisition, high resolution three dimensional experiments lately redefined state art. Tackling spatial temporal requires high‐performance computing approach. We developed multi‐ graphical processing units (GPU) application resolve anisotropic elastodynamic that relies on conservative scheme simulate, few seconds, fields for domains involving more than 1.5 billion grid cells. present comprehensive analysis reducing number material parameters needed from ten four. Furthermore, emphasizes key governing poroelastic media. perform dispersion function dimensionless leading simple transparent relations. then benchmark our against an analytical plane solution. Finally, we several modeling experiments, including three‐dimensional simulation injection into medium. provide Matlab, symbolic Maple, GPU CUDA C routines reproduce main presented results. efficiency implementation makes it readily usable investigate high‐resolution scenarios practical applications.